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We have carried out a preliminary design and simulation of a single-electron resistive switch 
based on a system of two linear, paraUel, electrostatically-coupled molecules: one implementing 
a single-electron transistor and another serving as a single-electron trap. To verify our design, 
we have performed a theoretical analysis of this memristive device, based on a combination of 
ab-initio calculations of the electronic structures of the molecules and the general theory of single- 
electron tunneling in systems with discrete energy spectra. Our results show that such molecular 
assemblies, with a length below 10 nm and a footprint area of about 5 nm^, may combine sub-second 
switching times with multi-year retention times and high (> 10^) ON/OFF current ratios, at room 
temperature. Moreover, Monte Carlo simulations of self-assembled monolayers (SAM) based on such 
molecular assemblies have shown that such monolayers may also be used as resistive switches, with 
comparable characteristics and, in addition, be highly tolerant to defects and stray offset charges. 



I. INTRODUCTION 

Recently, a substantial progress was made in the fab- 
rication of two-terminal "memristive" devices (includ- 
ing bistable "resistive" or "latching" switches) based on 
metal oxide thin-films, whose bistability is apparently 
based on the reversible formation/dissolution of conduct- 
ing filaments — see, e.g., recent reviews [TH5]. However, 
scaling of resistive memories and hybrid CMOS/nano- 
crossbar integrated circuits |T], based on such switches, 
beyond the 10 nm frontier may still require more repro- 
ducible devices based on other physical principles. One 
possibility here is to use a molecular version of single- 
electron switches [4]. 

Such a switch, schematically shown in Fig. la, is 
a combination of two electrostatically-coupled devices: 
a "single-electron trap" and a "single-electron transis- 
tor" [S] placed in parallel between two electrodes. (It 
will be convenient for us to call these electrodes the 
"drain/control" and the "source" — see Fig. la.) When 
the charge state of the trap island is electroneutral (Q = 
0), the Coulomb blockade threshold Vc of the transis- 
tor is large (Fig. lb), so that at all applied voltages with 
\V\ < Vc the transistor carries virtually no current — the 
so-called OFF state of the switch. As soon as the volt- 
age exceeds a certain threshold value < Vc, the rate 
of tunneling into the single-electron trap island increases 
sharply (Fig. Ic), and an additional elementary charge q 
(either a hole or an electron) enters the trap island from 
the source electrode, charging it to Q — q = ±e. The 
electrostatic field of this charge shifts the background 
electrostatic potential of transistor's island and as a re- 
sult reduces the Coulomb blockade threshold of the tran- 
sistor to a lower value VA. This is the ON state of 
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the switch, with a substantial average current flowing 
through the transistor a.t V > Vq. The device may be 
switched back into the OFF state by applying a reverse 
voltage in excess of the trap-discharging threshold \V^\. 
As was experimentally demonstrated for metallic, low- 
temperature prototypes of the single-electron switch [5], 
its retention time may be very long [7]. However, for that 
the scale e^/2C of the single-electron charging energy of 
the trap island, with effective capacitance C, has to be 
much higher than the scale of thermal fluctuations, k^T. 
For room temperature, this means the need for few-nm- 
sized islands [51, and so far the only way of reproducible 
fabrication of features so small has been the chemical 
synthesis of suitable molecules — see, e.g., [HIS]. 

Transport properties of single molecules, captured be- 
tween two metallic electrodes, have been repeatedly stud- 
ied by several research groups for more than a decade — 
see, e.g., [TOHTB] . The practical use of such devices in 
VLSI circuits is still impeded by the unacceptably low 
yield of their fabrication and large device-to-device vari- 
ability. The main reason of this problem is apparently 
the lack of atomic-scale control of the contacts between 
the molecules and the metallic electrodes. In addition, in 
three-terminal single-molecule devices that can work as 
transistors [13l [M] , there is an additional huge challenge 
of reproducible patterning of three-electrode geometries 
with the required sub-nm precision. These two chal- 
lenges make single-molecule three-terminal stand-alone 
devices rather unlikely candidates for post-CMOS VLSI 
circuit technology. However, we believe that for resis- 
tive memories and CMOS/nano hybrids based on nano- 
crossbars (Fig. Id pTj), these challenges may be met. 
Indeed, such crossbars use two-terminal crosspoint de- 
vices, so that their only critical dimension (the distance 
between the two electrodes) may be precisely controlled 
by layer thickness. In addition, if such devices are based 
on self-assembled monolayers (SAMs), the large num- 
ber N of molecules in a single device may mitigate the 
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negative effects of interfacial and other uncertainties of 
single molecules [II]. (Since the electrode footprint of 
a quasi-linear functional molecule stretched between the 
electrodes may be very small, N may be as high as ^ 10^ 
even for sub-lO-nm-scale devices.) 

The goal of this paper is to describe the results of the 
design and ab-initio calculations of basic properties of 
a molecular resistive switch and SAMs based on such 
molecules. Our basic design is described, and its physics 
is discussed in the next section (Sec. II). In Sec. Ill we 
formulate a theoretical model which allows an approx- 
imate but reasonably accurate numerical simulation of 
electron transport properties of this device. The results 
of simulation of our most promising switch version are 
described in Sec. IV. In Sec. V we describe our approach 
to simulation of SAM layers consisting of resistive switch 
assemblies, and the results of these simulations. Finally, 
in Sec. VI we summarize our results and discuss the nec- 
essary further work towards the practical implementation 
of reproducible resistive switches. 
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FIG. 1. (a) The traditional version of the single-electron re- 
sistive switch, (b) its 7 — curve (schematically) , and (c) the 
ON/OFF switching rates of the device, calculated using the 
orthodox theory of single-electron tunneling [TH], for /C = 
20k^T. The inset in panel (c) schematically shows the switch- 
ing between charge states of the trap resulting from repeated 
voltage sweeps with a rate Fo = \dV/dt\ /(e/C) ^ F,-. (d) 
Nano-crossbar with resistive switches as crosspoint devices. 



II. RESISTIVE SWITCH DESIGN 

Our initial design [20] of the molecular resistive 
switch was based on oligophenyleneethynylene (OPE) 



chains as tunnel barriers and diimide (namely, pyromel- 
litdiimide, naphthalenediimide, and perylenediimide) 
groups as trap and transistor islands. However, already 
the first quantitative simulations have shown that the 
relatively narrow HOMO-LUMO bandgap of the OPE 
chains (of the order of 1.5 eV [21 ) cannot provide a tun- 
nel barrier high enough to ensure sufficiently long elec- 
tron retention times in traps with acceptable lengths. As 
a result, we have concluded that alkane chains (CH2- 
CH2-...), with a bandgap of ~ 9 eV [22, are much bet- 
ter candidates for resistive switches. There is also sub- 
stantial experience in the chemical synthesis of molecular 
electronic devices and SAMs using such chains as tunnel 
barriers [HIS]. 

Figure 2 shows a possible realization of such a device, 
based on benzene-benzobisoxazole and naphthalenedi- 
imide acceptor groups (playing the role of single-electron 
islands), and alkane chains. In contrast with the usual 
( "orthodox" ) design of the trap j4, ^ i23j , where a long 
charge retention time is achieved by incorporation of sev- 
eral additional single-electron islands into the trap charg- 
ing path (Fig. la), the long alkane chain used in the 
molecular trap has a band structure which enables its 
use simultaneously in two roles: as a tunnel barrier as 
well as a replacement of intermediate islands. 
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FIG. 2. Our final version of a molecular resistive switch fea- 
turing an alkane-naphthalenediimide single-electron trap elec- 
trostatically coupled to an alkane-benzobisoxazole-benzene 
single-electron transistor. 

In order to explain this novel approach, let us first re- 
view the role of intermediate islands in the conventional 
design of the trap (Fig. la). If a single-electron island is 
so large that the electron motion quantization inside it is 
negligible, its energy spectrum, at a fixed net charge Q, 
may be treated as a continuum. Elementary charging of 
the island with either an additional electron or an addi- 
tional hole raises all energies in the spectrum by e^/2C, 
where C is the effective capacitance of the island [5J [33] ■ 
As a result, the continua of the effective single-particle 
energies of the system for electrons and holes are sepa- 
rated by an effective energy gap e^/C — essentially, the 
"Coulomb gap" Pl]. If this gap is much larger than /cbT 
at applied voltages V close to the "energy-equilibrating" 
voltage Ve (see the middle panel of Fig. 3a), it may 
ensure a very low rate F^ of single-charge tunneling in 
either direction and hence a sufficiently long retention 
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time tr = 1/Tr of the trap. The energy gap may be sup- 
pressed by applying sufficiently high voltages ~ e/C of 
the proper polarity, enabling fast switching of the device 
into the counterpart charge state — see the left and right 
panels of Fig. 3a, and also Fig. Ic. 
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FIG. 3. Schematic single-particle energy diagrams of (a) the 
usual single-electron trap shown in Fig. la (for the sake 
of simplicity, with just one intermediate island) and (b) the 
molecular trap analyzed in this work (Fig. 2), each for three 
values of the applied voltage V. Occupied/unoccupied energy 
levels are shown in black/green. (The dotted green/black line 
denotes the energy level of the "working" orbital that is either 
empty or occupied during the device operation, defining its 
ON/OFF state.) Horizontal arrows show (elastic) tunneling 
transitions, while vertical arrows indicate inelastic relaxation 
transitions within an island, a molecule, or an electrode. 

In the molecular single-electron trap shown in Fig. 2, 
the "energy-equilibrating" voltage Vg aligns the Fermi en- 
ergy of the source electrode with the lowest unoccupied 
level of the acceptor group that is, by design, located 
in the middle of the HOMO-LUMO gap of the alkane 
chain. As a result, an electron from the source electrode 
may elastically tunnel into the group only with a very low 
rate F^ — see the middle panel in Fig. 3b. The reciprocal 
process (at the same voltage) may be viewed as electron 
tunneling from the highest occupied molecular orbital of 
the singly-negatively charged molecule. (To simplify the 
terminology, in the reminder of the paper we call this 
molecular orbital the "working orbital" (indexed W) , in- 
stead of HOMO or LUMO of the counterpart molecular 
ions, to make the name independent of the charge state 
of the device.) The energy-balance condition of both pro- 
cesses is similar, and may be expressed via the effective 
single-particle energy ew of the working orbital |25) : 



where Egr{n) is the ground-state energy of the molecu- 
lar ion with n electrons. (In the case of singly-negative 
ion we are discussing, n = uq + 1 |26| . where hq is the 
number of protons in the molecule.) In this notation, the 
energy-balance (level-alignment) condition, which defines 
the voltage Ve, is 



ew = W - ejVe, 



(2) 



ew = ^E{n) = Egr{n) - Egr{n - 1), 



(1) 



where W is the workfunction of the source electrode ma- 
terial, and 7 is a constant factor imposed by the geometry 
of the junction; < 7 < 1. (Its physical meaning is the 
fraction of the applied voltage, which drops between the 
trapping island and the source electrode.) At the charg- 
ing threshold voltage V^, energy ew becomes aligned 
with the valence band edge of the chain, allowing for a 
fast charging of the molecule — see the left panel in Fig. 
3b. Similarly, as shown on the right panel in Fig. 3b, 
at VLj. this energy becomes aligned with the conduction 
band edge of the chain, allowing for a fast discharging of 
the molecule. 

As an example. Fig. 4a shows the atomic self- 
interaction corrected (ASIC) [57] Kohn-Sham elec- 
tron eigenenergy spectrum ef^^''^{n) of the alkane- 
naphthalenediimide molecule used as our final trap de- 
sign (Fig. 2), with the net charge Q{n) — —e{n — uq) — 
— e, as a function of the applied voltage V. (Here i is 
the spin-orbital index; see Sec. Ill below.) Point col- 
ors in Fig. 4a crudely represent the spatial localization 
of the orbitals: blue corresponds to their localization at 
the trapping (acceptor) group, while red marks the lo- 
calization at the alkane chain's part close to the source 
electrode. Figure 4b shows the probability density of the 
working orbital 4'w^'^ ~ ^no^+i('^o + 1) of the molecu- 
lar trap, integrated over the directions perpendicular to 
the molecule's axis, with blue lines corresponding to the 
probability density at the most negative applied voltage. 
At the equilibrating voltage Ve ~ 2.2 V, the working 
orbital is well localized at the acceptor group, and is iso- 
lated from the source electrode by a ^ 4.5-eV-high energy 
barrier. However, as Fig. 4b shows, the probability den- 
sity of the orbital decays into the alkane group rather 
slowly, with the exponent coefficient (3 w OAa^^, corre- 
sponding (in the parabolic approximation of the disper- 
sion relation) to the effective electron mass nief ~ O.lmo 
P5] . As a result, a long (~ 5 nm) alkane chain is needed 
to ensure an acceptable retention time of the trap. (The 
2-nm free-space separation between the other side of the 
molecule and the control/drain electrode, shown in Fig. 
2, is quite sufficient for preventing electron escape in that 
direction.) 

At a sufficiently high forward/reverse bias voltage, 
the working orbital energy ew crosses into the conduc- 
tion/valence band of the alkane chain, so that the or- 
bital partly hybridizes with the states localized near the 
source electrode interface, described by the rise of IV'wP 
at larger values of z — see Fig. ^p. This rise enables 
fast electron tunneling to/from the source electrode, i.e. 
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FIG. 4. ASIC density-functional-theory (DFT) results (cor- 
rected for level "freezing" at high positive and negative volt- 
ages — see Sec. Ill for details) for the singly-negatively 
charged alkane-naphthalenediimide trap molecule, (a) Kohn- 
Sham energy spectrum as a function of the applied voltage V, 
with colors representing the spatial localization of the corre- 
sponding orbitals — see the legend bar on the right. The ver- 
tical lines mark voltage values !/<_, Ve, and corresponding 
to the left, middle, and right panels of Fig. 3b. The dashed 
lines labeled E-p{s) and E-p{c) show the Fermi levels of the 
source and control/drain electrodes whose workfunction was 
assumed to equal 5 eV. (b) The "working" orbital's proba- 
bility density, integrated over the directions perpendicular to 
molecule's axis, for a series of applied voltages V — see the 
legend bar on the right of the panel. 



a fast switching of the device to the counterpart charge 
state, in a manner similar to that of the conventional 
single-electron trap, as shown schematically on the left 
and right panels of Fig. [3]3. Thus the long molecular 
chain, with a sufficiently large HOMO-LUMO gap, may 
indeed play the roles of both the tunnel junction and in- 
termediate islands of the "orthodox" single-electron trap. 



For the design of the second component of the switch, 
the molecular single-electron transistor, the most impor- 
tant challenge is to satisfy the ON and OFF state current 
requirements. In particular, the ON current should not 
be too large to keep the power dissipation in the circuit 
at a manageable level, but simultaneously not too small, 
so that the device output does not vanish in the noise of 
the sense amplifier (for memory applications [29, 30 ) or 
the CMOS invertor (in hybrid logic circuits |^). Also, 
the ON/OFF current ratio should be sufficiently high to 
suppress current "sneak paths" in large crossbar arrays 
[nnilSl]- In addition, the transistor molecule should be 
geometrically and chemically compatible with the trap 
molecule, enabling their chemical assembly as a uni- 
molecular device, with their single-electron island groups 
sufficiently close to provide substantial electrostatic cou- 
pling. (Without it, the charge of the trap would not 
provide a substantial modulation of the transistor cur- 
rent.) At the same time, the molecules must not be too 
close, in order to prevent a parasitic discharge of the trap 
via electron cotunneling through the transistor into one 
of the electrodes. The chemical compatibility strongly 
favors the use of similar chains as the transistor's tunnel 
junctions. 



We have analyzed several alkane-chain based tran- 
sistor devices with naphthalenediimide, perylenediimide 
and benzobisoxazole acceptor groups as transistor is- 
lands. However, in all these cases the long alkane chains, 
needed to match the lengths of the transistor and trap 
molecules, make ON currents too low. Finally, we have 
decided to use an unusually long (~ 4.3-nm) benzene- 
benzobisoxazole |33j island group — see Fig. [2] and Fig. 
[5]d. Figure [5^ shows the Kohn-Sham electron eigenen- 



ergy spectrum e 



ASIC 



-1) of this molecule as a function 
of voltage V (where Iq is the number of protons in the 
transistor molecule). Blue/red colored points correspond 
to the orbitals localized at the left/right alkane chain, 
while green color points denote the orbitals extended over 
the whole acceptor group. This extension is clearly vis- 
ible in Fig. 5b, which shows the probability density of 
the working orbital e 



ASIC _ ^ASIC 

w 



£(0+1(^0 + 1) of the transis- 
tor molecule. During transistor operation, the tunneling 
electron may populate any of several group-localized or- 
bitals, resembling the operation of the usual (metallic) 
single-electron transistor. As a result of such island ex- 
tension, alkane chains of the transistor could be substan- 
tially shortened, to ~ 1.5-nm-long C11H25, enabling low 
but still acceptable ON currents of the order of 0.1 pA, 
even if a small (0.25-nm) vacuum gap between the alkane 
chain and the source electrode is included into calcula- 
tions to give a phenomenological description of the ex- 
perimentally observed current reduction due to unknown 
interfacial chemistry |25j . 
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FIG. 5. ASIC results for the single-negatively charged 
benzene-benzobisoxazole transistor molecule, (a) Kohn-Sham 
energy spectrum as a function of the applied voltage V, with 
colors representing the spatial localization (within the junc- 
tion) of the corresponding orbitals — see the legend bar on 
the right. The dashed lines labelled E-pi^s) and Ep{c) show 
the Fermi levels of the source and control/drain electrodes 
whose workfunction was assumed to equal 5 eV. (b) Proba- 
bility density of the working orbital, integrated over the di- 
rections perpendicular to the molecular axis, for a series of 
applied voltages — see the legend bar on the right. 



III. THEORETICAL MODEL AND 
APPROXIMATIONS 

Each molecule used in our device has a discrete set 
of possible excited states, and hence the electron trans- 
port is not limited to a single channel. In order to 
take into account all of these channels, we have used 
the "quasi-single-particle approximation" whose simplest 
version had been first formulated by Averin and Korotkov 
for semiconductor quantum dots |18l 134) and which was 
recently generalized [25j to be more applicable to molec- 



ular structures. In this approximation, the energy of an 
arbitrary state k = {n, i} of the molecule equals 

Ek = Egr{n) +^ei{n)p, -^e^{72){l - Pi), (3) 



where coefficients ei{n) have the physical meaning of 
single-particle excitation energies of an n-electron ion, 
and numbers pi (equal to either or 1) are the single- 
particle energy level occupancies. The condition of elastic 
tunneling, leading to a transition between states k and 
k' , is given by the natural generalization of Eqs. (fll) and 



£k^k' = W - ejV, 



(4) 



where the single-electron recharging/excitation energy is 
now defined as 



£k^k' 



Ek — Ehi ■ 



(5) 



Because of the large size and complexity of the 
molecules used in our design, the only practical way to 
calculate their electronic structure is to use a software 
package (such as SIESTA [3S] ) , based on the density- 
functional-theory (DFT) [37], which may provide a rea- 
sonably accurate ground state energy E^^'^{n) and a 
single-particle spectrum eP^""" (n) at a fraction of the com- 
putational cost of more correct ab-initio methods. Unfor- 
tunately, for such strongly correlated electronic systems 
as molecules considered in this paper, results obtained 
using standard DFT software packages [35] have signifi- 
cant self-interaction errors [39j . 

We believe the source of such errors is that the ap- 
proximate treatment of the exchange-correlation term in 
the Kohn-Sham Hamiltonian does not completely can- 
cel the self-interaction energy present in the "Hartree 
term" of the Hamiltonian [40 . Indeed, the standard 
DFT approach leads to errors, in the key energies ^ 
and ([5|, of the order of the single-electron charging 
energy e^/2C, where C is the effective capacitance of 
the island group — see Appendix A for details. This 
error may be rather substantial; for example in the 
naphthalenediimide-based trap molecule (Fig. 2), it is 
approximately equal to 1.8 eV. For this reason, the elec- 
tron affinity Egr{no + 1) — Egr{nQ), calculated using the 
LSDA DFT for the singly-negatively charged ion of the 
molecular trap, is significantly (by ~ 3.2 eV) larger than 
the experimental value of similar molecules ,41, .42 . The 
LSDA energies may be readily corrected to yield a much 
better agreement with experiments (see Table 1 in Ap- 
pendix A), however, it is not quite clear how such a the- 
ory may be used for a self-consistent calculation of the 
corresponding working orbital ■(/'i^(r). 

We have found that a significant improvement may 
be achieved by using the recently proposed Atomic Self- 
Interaction Corrected DFT scheme (dubbed ASIC [27]) 
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implemented in a custom version of the SIESTA software 
package. For the molecules that we have considered here, 
this approach gives the Kohn-Sham energy e„o+i(?io + 1) 
very close to the experimental electron affinity. However, 
we have found that using even this advanced approach for 
our task faces two challenges. 

First, the algorithm gives (at least for our molecular 
trap states with n = uq + 1 and n = Uq + 2 electrons) 
substantial deviations from the relation ew = £n+i{n+l) 
for n = riQ (which has to be satisfied in any exact theory 
[571 03]), with the ground energy difference (Q^ close to 
the LSDA DFT results. This means that Eq.^ cannot 
be directly used with the ASIC results; instead, for the 
electron transfer energy between adjacent ions n and n— 1 
we have used the following expression: 



-ASIC 



(n). 



(6) 



This relation implies that the differences £^^^'^( 



-ASIC 



(n) describe all possible single-particle excitations 
within the acceptor group, if the index i' is restricted 
to orbitals localized on the group. (Other orbitals, local- 
ized on the alkane chain are irrelevant for our calculations 
since they do not contribute to the elastic tunneling be- 
tween the molecular group and the electrode.) 

In order to appreciate the second problem, look at Fig. 
6 which shows the voltage-dependent Kohn-Sham spec- 
tra of the singly-negatively charged molecular trap, cal- 
culated using the ASIC SIESTA package for T > K. 
Notice that above voltage Vt ~ 13V, and below volt- 
age Vf « — 7V, the eigenenergy spectrum is virtually 
"frozen". (The LSDA SIESTA gives similar results.). 
As explained in Appendix B using a simple but rea- 
sonable model (similar to that used in Appendix A), at 
V > Vt such "freezing" originates from the spurious self- 
interaction of an electron whose wavefunction cloud is 
gradually shifted from the top occupied orbital of the va- 
lence band of the chain, with energy Sy , into the initially 
empty group-localized orbital with energy Sw+i- (A sim- 
ilar freeze at voltages V < V/, is due to the spurious 
gradual transfer of the electron wavefunction cloud from 
the working orbital, localized at the acceptor group, with 
energy ew, to the lowest orbital of the conduction band 
of the chain, with energy Sc-) It is somewhat surprising 
that this spurious effect (which should not be present in 
any consistent quantum-mechanical approach — see Ap- 
pendix B) is so strongly expressed in the ASIC version 
of the SIESTA code, which was purposely designed to 
get rid of the self-interaction in the first place. Being no 
SIESTA experts, we may only speculate that the nature 
of this artifact is related to the smoothing of the deriva- 
tive discontinuity present in the ASIC method as the elec- 
tron number passes through an integer value, which is 
mentioned in [57| — see also Fig. 7 in that paper. 

Fortunately, there is a way to correct this error very 
substantially by following the iterative process of self- 
consistent energy minimization within ASIC SIESTA. In- 
deed, for a fixed temperature T > K (when the program 
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FIG. 6. The Kohn-Sham spectra of the singly-negatively 
charged molecular trap, calculated with the ASIC SIESTA 
at r = 10 K. At voltages below V/ ~ —7 V, the spectrum is 
virtually frozen due to a spurious gradual shift of the highest- 
energy electron from the "working" orbital (with energy ew, 
shown with a solid blue line) localized on the acceptor group, 
to the lowest orbital (with energy Ec, shown with a solid red 
line) of the conduction band of the alkane chain. As a re- 
sult, the calculated spectrum is virtually voltage-insensitive 
("frozen"). In the voltage range Vt < V < Vt ^ 13 V, ASIC 
SIESTA gives apparently correct solutions, with the work- 
ing orbital ew fully occupied, and the next group-localized 
orbital (with energy ew+i, the dashed blue line) completely 
unoccupied. However, a.t V > Vt the package describes a 
similar spurious gradual shift of the highest-energy electron 
from the highest level of the valence band of the chain to 
orbital ew+i, resulting in a similar spectrum "freeze". The 
spectrum evolution, calculated after the (approximate) cor- 
rection of this spurious "freezing" effect, is shown in Fig. 4a 
above. 



automatically populates molecular orbitals in accordance 
with the single-particle Fermi-Dirac statistics) and volt- 
ages y > 14 w 13 V and 1/ < F/ « -7 V, its iterative 
process converges to a wrong solution with the energy 
levels frozen at their Vt and V/ values, as is discussed 
above — see Fig. 6. However, if the temperature in that 
program is fixed at T = K, its iterative process ends 
up in quasi-periodic oscillations between different solu- 
tions — most of them with frozen levels (just like in Fig. 
6), but some of them with the group localized energies 
like the working orbital energies ew, £w+i and the va- 
lence/conduction band edge energies e^, ec close to their 
expected (unfrozen) values. (Those values were obtained 
by a linear extrapolation of their voltage behavior calcu- 
lated at V/ <V < Vt.) Since such a solution is repeated 
almost exactly at each iterative cycle (see the vertical 
boxes in Fig. 7), we believe that it is close to the cor- 
rect solution expected from the self-consistent quantum- 
mechanical theory — see Appendix B. These approxi- 
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mate solutions were used in our calculations both above 
Vt and below V/; they are illustrated in Fig. 4a, where 
we have substituted the incorrect "frozen" solutions for 
T > K with solutions for T = K, with ew ~ 
ew+i ~ £w+i s-i^d Sc ~ at y < V/ or « e^* at 
V > Vt- Let us emphasize that the approximate nature 
of these solutions may have affected our calculations (we 
believe, rather insignificantly), only at F > Vt w 13 V 
and V < V{ ~ —7 V, i.e. only the device recharging 
time results, but not the most important retention time 
calculations for smaller voltages — see Fig. 9 below. 
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FIG. 7. The Kohn-Sham energy spectrum of our trap 
molecule, as calculated by successive iterations within ASIC 
SIESTA for r = K and 1/ = 14.9 V, i.e. above the threshold 
voltage Vt ~ 13 V. Vertical boxes mark the apparently cor- 
rect solutions with energies of the working orbital (ew^^^), 
the next group-localized orbital {eyv+i)^ ^nd the highest or- 
bital of the valence band of the alkane chain (e^^^*"^) all close 
to their respective values e^, £w+i a-nd ej}* obtained by a 
linear extrapolation of their voltage dependence calculated at 
Vt <V <Vt- Just like in Figs. 4a, 5a and 6, point colors rep- 
resent the spatial localization of the corresponding orbitals. 
Lines are only guides for the eye. 

With the electron orbitals and eigenenergies calcu- 
lated, we have described dynamics of both the trap and 
the transistor, just as in our first work [3S], by a set 
of master equations for state probabilities |18j, which 
are valid because of the incoherent character of single- 
electron tunneling to/from the continuum of electronic 
states in metallic electrodes [S]. Moreover, for the inelas- 
tic relaxation rates Fjnei and the rates F<_ and F_j. of the 
elastic tunneling between the molecular group and the 
metallic electrodes (see arrows in Fig. 3), the following 
strong inequality. 



incl 



(7) 



is well fulfilled. (Indeed, the rates Finei are crudely of the 
order of 10^^ 1/s in both molecules and metals. On the 



other hand, our results, described in Sec. IV below, yield 
transistor currents I ^ lO""'^'^ A, meaning that F<_ and 
F^ are of the order of I /e ^ 10^ 1/s in the transistor; 
the rates are even much lower than that in the trap — see 
Fig. 9b.) Relation Q allows us to account only for the 
tunneling events starting from thermal equilibrium, and 
ensures that the rates Fjnei drop out of the calculations. 

In comparison with |25j , one more new element of this 
work is the electrostatic coupling between the trap and 
the transistor which features similar but much more fre- 
quent single-charge transitions. This rate hierarchy al- 
lows the trap to be described by averaging rates F of 
tunneling events in it over a time interval much longer 
than the average time period between tunneling events 
in the transistor, but still much shorter than 1/F. These 
average rates may be calculated as 



(8) 



for electron tunneling from the source into the trap 
molecule, and 



(9) 



for the reciprocal event. Here (t„(Z) are the conditional 
probabilities of certain charge states / of the transistor 
island provided that the trap is in the n-electron charge 
state (with n equal to either rip or no + 1), while wi^±{n) 
are the total rates of single electron tunneling between 
the trap in its initial charge state n and the source elec- 
trode. These rates have been calculated using Eq. (11) in 
|25) , with an extra index I added to account for the tran- 
sistor's state. The conditional probabilities cr„(/) satisfy 
the usual normalization condition: 



(10) 



and (together with the dc current / fiowing through the 
transistor) have been calculated as in }25| . by combining 
the master equations of single-electronics [151 [M] with 
ab-initio calculations of molecular orbitals and spectra, 
and the Bardeen formula [4^ for tunneling rates. 

The electrostatic interaction between the two 
molecules is taken into account by an iterative incorpo- 
ration of the numerically calculated Coulomb potential 
created by both molecules (as well as by the series of 
their charge images in the metallic electrodes of the 
system, which we have assumed to be plane, infinite 
surfaces — see Fig. [s]) into the Kohn-Sham potentials. 
From the elementary electrostatics, this potential may 
be expressed as 
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FIG. 8. A schematic view of charge densities participating in 
Eq. ([11}. 
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where Ps(ro), Pc(ro) are the total charge distributions 
of the molecule under analysis and the complementary 
molecule, and ps(rj), Pc(jj) are the corresponding charge 
images in the source (j > 0) and the drain {j < 0) elec- 
trodes — see Fig. [Sj The first term in Eq. (11) is the 



potential created by the complement molecule, the sec- 
ond term describes the potential of the infinite set of 
charge images of both molecules in the source and con- 
trol/drain electrodes, and the third term is the potential 
created by the applied source-drain voltage. At the 0-th 
iteration, the first two terms are taken equal to zero. For- 
tunately, the iterations give rapidly converging results, so 
that there was actually no need to go beyond the second 
iteration — see Fig. [9j 

Another important change introduced into our calcu- 
lations is that the charge transfer rates (see Fig. 4c in 
|25j ) are calculated in a simpler way. Namely, one of the 
key conditions of validity of the Bardeen formula for the 
tunneling matrix elements 



2m 



din 

'' dz 



dz 



dS 



(12) 



(where ipi is the molecular orbital, il)s,c are the wave- 
functions of electrons located inside the source or con- 
trol/drain electrodes, and S is an arbitrary surface sepa- 
rating the single-electron island from the corresponding 
electrode) is that the result given by Eq. ( 12 ) is inde- 
pendent of the position of the surface S. Due to elec- 
trostatic screening of the electric field by the electrode, 
the Kohn-Sham potential becomes very close to the vac- 
uum potential at just a few Bohr radii ob away from the 
molecule's last atom. Therefore, if the surface S is se- 
lected inside the vacuum gap between the molecule's end 




1 
V (Volts) 

FIG. 9. Differences between the energies of the working or- 
bital of the molecular trap, calculated with ASIC SIESTA at 
the fc-th and {k — l)-th iterations of the Coulomb interac- 
tion potential, given by Eq. (Ill, as functions of the applied 
voltage. 



and the electrode surface, the effect of the molecule on 
wavefunctions 'ips,c is negligible with good accuracy (cor- 
responding to a fraction of one order of magnitude in the 
resulting current). Hence, these wavefunctions may be 
calculated analytically to describe the usual exponential 
ID decay into vacuum, instead of a numerical solution of 
a Schrodinger equation, as it had been done in |25j. 

On top of the ab-initio calculation scheme, we have 
performed the following check of the component molecule 
spacing. As was mentioned in Sec. II, the molecules have 
to be placed sufhciently far from each other to prevent a 
parasitic discharge of the trap via the elastic cotunncling 
through the transistor island, into one of the electrodes. 
This effect may be estimated using the following formula 



pcot 



2A2 + 4e2 • 



(13) 



where A is the matrix element of electron tunneling be- 
tween the trap and the transistor islands (that exponen- 
tially depends on the distance dp between the molecules, 
see Fig. [2]), e is the difference between eigenenergies of 
these states, and is the rate of tunneling between the 
transistor island and its electrodes. We have found that 
in order for the cotunneling rate F™* to be below the 
retention rate of the trap F^, the distance dp should not 
be lower than ~ 1.5 nm. Such relatively large separa- 
tion justifies separate DFT calculations of the electronic 
structures of the trap and the transistor (with the molec- 
ular geometry of each component device initially relaxed, 
using LSDA SIESTA gT]), related only by their electro- 
static interaction described by Eq. (11). 



9 



IV. SIMULATION RESULTS FOR A SINGLE 
RESISTIVE SWITCH 



Figure [TO] shows our main results for the resistive 
switch shown in Fig. [2] for temperature T — 300 K. 
They include the dc I — V curves of the transistor, plot- 
ted in Fig. 10 i for both charge states of the trap, and 



10 1), and the ON current / 



0.2 



the rates of transitions between the neutral and single- 
negatively charged states of the trap, with and without 
the account of the transistor effect on the trap molecule 
(Fig. [lOjs). The plots show that the resulting I—V curves 
fit our initial specifications rather well, with a broad volt- 
age window (from ~ 2.0 to ~ 2.5 V) for the trap state 
readout, a large ON/OFF current ratio within that win- 
dow (inset in Fig 
pA. 

Figure [10|d shows that the trap features a high reten- 
tion time, Tr > lO^s, for both charge states, within a 
broad voltage range, — 2V < V < -1-5 V. (It is some- 
what surprising how little is the trap retention affected 
by the electrostatic "shot noise" generated by fast, quasi- 
periodic charging and discharging of the transistor island, 
which is taken into account by our theory.) The range in- 
cludes point V = 0, so that the device may be considered 
a nonvolatile memory cell. 

At the same time, the device may be switched between 
its states relatively quickly by applied voltages outside of 
this window. The price being paid for using alkane chains 
with their large HOMO-LUMO gap is that the voltages 
necessary for fast switching are large — they must align 
the valence or conduction band of the alkane chain with 
the group-localized working orbital — see Fig. [4^. 



V. SAMS OF RESISTIVE SWITCHES 

Probably the largest problem of molecular electron- 
ics [8, 9 is the low reproducibility of interfaces between 
molecules and metallic electrodes. However, recent re- 
sults [IH] indicate that this challenge may be met at least 
for self-assembled monolayers (SAMs) encapsulated us- 
ing special organic counter-electrodes. This is why we 
have explored properties of SAMs consisting of square 
arrays of N x N resistive switches described above — see 
Fig. [TT] In order to increase the tolerance of the resulting 
SAM devices to self-assembly defects and charged impu- 
rities, it is beneficial to place the component molecular 
assemblies (Fig. 2) as close to each other as possible, say 
at distances comparable to that (^1.5 nm) between the 
trap and transistor. In this case, the Coulomb interac- 
tions between the component molecules are very substan- 
tial, and properties of the system have to be calculated 
taking these interactions into account. 

A system oi N x N resistive switches has 2A^ x 2A^ 
single-electron islands and hence at least 2^^^^^ possi- 
ble charge states, which would require solving that many 
master equations for their exact description. Even for 
relatively small N, this approach is impracticable, and 




5 
V (Volts) 

FIG. 10. (a) Calculated dc I — V curves of the transistor for 
two possible charge states of the trap molecule. The inset 
shows the ON/OFF current ratio of the transistor on a semi- 
log scale, within the most important voltage interval, (b) The 
trap switching rates, calculated with (solid lines) and without 
(dashed lines) taking into account the transistor's back action, 
as functions of the applied voltage. Red dashed lines on panel 
(b) show the trap switching rates calculated without the level 
"freezing" correction (Sec. Ill and Appendix B) and without 
taking into the account the transistors' back action. 



virtually the only way to explore the properties of the 
system is to perform its Monte Carlo simulations O [35] . 
In this method a random number generator is used twice 
for each state change: first, to calculate the random time 
of some state change (which obeys Poissonian statistics) , 
and second, to calculate the charge transition type (if 
several transitions are possible simultaneously). The pro- 
cedure requires a prior calculation of rates of transitions 
between all pairs of charge configurations which differ by 
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FIG. 11. Schematic view of a 5 x 5-switch SAM sandwiched 
between two electrodes. 



one single-electron tunneling event. 
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FIG. 12. Effect of a single charge of a trap molecule on the 
electron affinity of another molecule, located at distance r 
without and with the account of the electric field screening 
by the common metallic electrodes. 

As was discussed above, the peculiarity of our partic- 
ular system is that it features two very different time 
scales: the first one (for our devices, Tt ~ e/IoN ^ 
10~^— 10~^ s) characterizes fast charge tunneling through 
single-electron transistors, and the second one corre- 
sponds to the lifetimes of trap states (r^ = l/T ^ 
10* — 10"'^ s). In order to gather reasonable statistics 
of the switching rates, our data accumulation time, for 
each parameter set, corresponded to the physical times 
of up to 10 s, i.e., included up to a million transition 
tunneling events in the system's transistors. 

As a check of the validity of the procedure, the Monte 
Carlo algorithm was first applied to a single resistive 
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FIG. 13. Trap tunnel rates as functions of the applied voltage 
for two quasi-similar nearest-neighbor charge configurations 
shown in the insets. 



switch, and it indeed gave virtually the same result as 
the master equation solution. We then used the ap- 
proach for a direct simulation of SAM fragments with 
two and more coupled resistive switches. As the frag- 
ment is increased beyond a 2 x 2 switch array, even the 
Monte Carlo method runs into computer limitations, be- 
cause of the exponentially growing number of the possible 
charge configurations. The calculations may be very sig- 
nificantly sped up by using the approximation in which 
each molecule's state affects the potential of only its near- 
est neighbors. This approximation has turned out to be 
very reasonable (Fig. 



12) and may be justified by the 
fact that metallic electrodes of the system substantially 
screen the Coulomb potential of the charges of distant 
molecules: the distance between the acceptor group cen- 
ters and the electrodes, d/2 « 4 nm, is of the same order 
as the 3-nm distance between the molecule and its next- 
next neighbors. In this nearest-neighbor approximation, 
each molecule (a trap or a transistor) is still affected by 8 
other molecules. To limit the number of the charge con- 
figurations even further, we have treated all "essentially 
similar" of them (having charge pairs at equal distances, 
irrespective of their angular position) as identical — see 
Fig. [13} 

Figure [M] shows the results of calculations, based on 
this approach, for a 5 x 5-switch SAM, of the total area 
close to 10 x 10 nm^. The switching and state read- 
out properties are very comparable with those of a sin- 



gle switch (Fig. 10), despite a significant mutual repul- 
sion between single electrons charging neighboring traps. 
In order to better understand why this repulsion does 
not have adverse effects on the operation of the SAM as 
a whole, we have calculated the correlation coefficients 
of charging of two molecules in the SAM as a function 
of the distance between them. At voltages above the 
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FIG. 14. Monte-Carlo simulated dc I — V curves of a 25- 
switch SAM. The top inset shows the fraction /3 of single- 
negatively charged traps, averaged over 40 sweeps of applied 
voltage between -8V and 13V. The bottom inset shows the 
ON/OFF current ratio averaged over the voltage sweeps, and 
its maximum sweep-to-sweep spread. 



transistor Coulomb blockade, transistor molecules switch 
their charge state fast and the correlation coefhcient K{r) 
between two transistor molecules may be calculated di- 
rectly from their time evolution records at a constant V . 
On the other hand, trap molecules have quasi-stationary 
charge states, so that the correlation between two trap 
molecules has to be calculated from a set of snapshots of 
their charge states (at some voltage of interest) taken at 
repeated, slow sweeping of the applied voltage through- 
out the whole voltage range. 



Figure 15 shows the resulting average correlation be- 
tween molecules (and its fluctuations) as a function of 
the distance between them in the 5 x 5-switch SAM. The 
charge states of neighboring traps are significantly an- 
ticorrelated, while the next-next neighbor charge states 
are positively correlated. This means that the switch- 
ing is due to a nearly-simultaneous entry of electrons 
into roughly every other trap [49 . This explains why 
in the top inset in Fig. [14] the average fraction of 
charged traps is close to 1/2. Thus the only adverse 
effect of the Coulomb interaction between individual re- 
sistive switches is the approximately two-fold reduction 



of the average ON current per device. Figure 16 presents 
a summary characterization of the SAM operation as a 
function of its size (and hence its area). 

The fact that even the fractional charging of traps in 
SAMs is sufficient for a very good modulation of their net 
current suggests that these devices should have a high tol- 
erance to defects and stray electric charges [5] . In order to 
verify this, we have carried out a preliminary evaluation 
of the defect tolerance by artificially fixing charge states 
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FIG. 15. The average correlation between two traps (green) 
and two transistors (blue) as a function of distance between 
them in a 5 X 5 device SAM. 
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FIG. 16. Summary of Monte Carlo simulations of SAMs of 
various area: (a) the average ON currents at voltages pro- 
viding certain ON/OFF ratios; (b) the average fraction of 
negatively charged traps at the equilibrating voltage Ve- 



of certain, randomly selected component molecules. The 
results, shown in Fig. [17] are rather encouraging, im- 
plying that the switches may provide the ON/OFF cur- 
rent ratios above 100 at defective switch fractions up to 
^ 10%, and at a comparable concentration of random 
offset charges. 
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number of random defects 

FIG. 17. Defect tolerance of the 5 x 5 SAM switch: ON 
current as a function of a number of molecules held artificially 
in a fixed, random charge state, at random locations, at the 
applied voltage values necessary to ensure a certain level of 
the ON/OFF current ratio. Error bars show the r.m.s. spread 
of results. 



VI. CONCLUSION 

Despite the problems with the description of single- 
electron charging in the density-functional theory, de- 
scribed in detail in Appendices A and B, we have man- 
aged to combine its advanced (ASIC) version to ana- 
lyze the possibility of using single-electron tunneling ef- 
fects in molecular assemblies for the implementation of 
bistable memristive devices ("resistive switches"). Our 
results indicate that chemically-plausible molecules and 
self-assembled monolayers of such molecules may indeed 
operate, at room temperature, as nonvolatile resistive 
switches which would combine multi-year retention times 
with sub-second switching times, and have ON/OFF cur- 
rent ratios in excess of lO'^. Moreover, we have obtained 
strong evidence that operation of the SAM version of the 
device is tolerant to a rather high concentration of de- 
fects and randomly located charged impurities. The ON 
current of a single device (~ 0.1 pA at F « 2 V) cor- 
responds to a very reasonable density (~ 4 W/cm^) of 
the power dissipated in an open SAM switch, potentially 
enabling 3D integration of hybrid CMOS/nano circuits 
[50] . (Note that the average power density in a cross- 
bar is at least 4 times lower because of the necessary 
crosspoint device spacing (Fig. [TJi); besides that, in all 
applications we are aware of, at least 50% of the switches 
(and frequently much more) are closed, decreasing the 
power even further.) 

However, even our best design (Fig. [2| still requires ad- 
ditional work. First, proper spatial positions of the func- 
tional molecules have to be enforced by some additional 



molecular support groups which have not been taken into 
account in our analysis yet. If the spacer groups fix- 
ing the relative spatial arrangement of the islands can 
be constructed from saturated molecular units similar to 
the alkane chains used to separate the islands from the 
electrodes, then the calculations presented here should 
be applicable to complete devices, but this expectation 
still has to be verified. 

Second, we feel that there is room for improvement 
in the choice of molecular chains used as tunnel barriers 
and intermediate islands. For example, the low calcu- 
lated effective mass, ruef ~ O.Itoq, of electrons tunneling 
along alkane chains makes it necessary to use rather long 
chains, despite their large HOMO-LUMO gaps (which, in 
turn, require large switching voltages — see Figs. [4{ [lO|. 
The use of a molecular chain with a higher rrief and a nar- 
rower gap would decrease switching voltages (and hence 
energy dissipation at switching), and also reduce the to- 
tal device length, resulting in shorter switching times (at 
the same charge retention). 

Third, the defect tolerance of SAM-based switches 
should be evaluated in more detail, for charged impu- 
rities located not only on the molecular acceptor groups, 
but also between them — say, inside the (still unspeci- 
fied) support groups. 

Finally, an experimental verification of our predictions 
looks imperative for the further progress of work towards 
practicable molecular resistive switches. 
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Appendix A: Single-electron charging correction 

Let us consider a simple but reasonable model of a 
well-conducting (say, metallic) island, of a size well above 
the Thomas-Fermi screening length, in which the single- 
electron addition energies are simply 

AEii) =K,~ e(/)„ (Al) 

where is i-th electron's kinetic energy (which, as well 
as the island capacitance C, is assumed to be independent 
of other electron state occupancies, but is an arbitrary 
function of i), and the second term describes the potential 
energy of that electron in the net electrostatic potential 
of all other charges, 
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bo^{t- 1) 



(A2) 



where (po is the background potential of the nuclei, and 
the second term is due to the previously added elec- 
trons. In this model the total ground-state energy of 
an n-electron ion (besides the electron-independent con- 
tributions) is 



^ K, - en4>o + - 1), (A3) 



so that the energy difference created by the last charging 
is 



AE{n) = Egr{n) - Egr{n - 1) 



= Kn- e(j)o + —{n - 1). 



(A4) 



On the other hand, in a hypothetical naive DFT 
theory, without the partial self-interaction corrections 
present in its LSDA, GGA and ASIC versions, the single- 
particle (Kohn-Sham) energies of ion n of this model are 
written as 



.DFT 



(n) ^ Ki- e(j)n, 0n = 00 - -^n 



(A5) 



For the calculation of the full ground-state energy of ion 
n, such generic DFT sums up these energies from i = I 
to j = n, adding the "double-counting correction" term 
[37] . in the Gaussian units equal to 



1 



(A6) 



where p(r) is the total electron charge density at point r. 
For our simple model, this correction is just e^ri^/2C, so 
that 



Comparing this result with Eq. (A4), we obtain the fol- 
lowing relation: 



AE{n) = AE°^^{n) 



e 

2C' 



(A9) 



Thus in the naive DFT theory, the single-electron ad- 
dition energy differs from the correct expression ( A4 ) by 



e /2C. Moreover, it does not satisfy the fundamental 



Eq. ([I]). Indeed, for i — n, Eq. (A9) gives the following 
result. 
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which, according to Eqs. ( A4 ) and ( A8 1 may be rewritten 
either as 



or as 
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This error is natural, because such DFT version ignores 
the fundamental physical fact that an electron does not 
interact with itself, even if it is quantum-mechanically 
spread over a finite volume. This difference can be- 
come quite substantial in small objects such as molecu- 
lar groups. For example. Table 1 shows the results using 
LSDA SIESTA calculations for two different ions of our 
trap molecule (Fig. [2]), with n = no 4- 1 and n = Uq + 2, 
where hq — 330 is the total number of protons in the 
molecule. The results show that the inconsistency de- 
scribed by Eq. (A12) is indeed very substantial and is 
independent (as it should be) of the applied voltage V 
in the range keeping the working orbital's energy inside 
the HOMO-LUMO gap of the alkane chain. The two last 
columns of the tables show the values of e^/2C, calcu- 
lated in two different ways: from the relation following 
from Eq. (|A5|): 
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and from the direct electrostatic expression 
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2C 



(A13) 



(A14) 



where (/)„ (r) is the part of the electrostatic potential, cre- 
ated by the electron of the n-th orbital of the n-th ion. 
The values are very close to each other and correspond 
to capacitance C « 4.5 x 10~^° F which a perfectly con- 
ducting sphere of diameter d « 0.8 nm would have. The 



14 



TABLE I. Columns 2 and 3: values of the single-electron 
transfer energy AE{n) for the trap molecule ions with n = 
no + 1 and n — no + 2 electrons, calculated in LSDA SIESTA 
and then self-interaction corrected as discussed in Appendix 
A, as functions of the applied voltage (Column 1). Columns 
3 and 4 list the values of parameter e^/C , calculated as dis- 
cussed in Appendix A. 



first empty 
group level 



-e{V- V, 



chain valence band 



Voltage V 
(V) 



AE{n) 
from Eq. 



(A9l (eV) 



AE{n) 
from Eq. 
(ATTI (eV) 



e72C 
from Eq. 
(AT3| (eV) 



e'/2C 
from Eq. 
(AT4| (eV) 



n = no + 1 


-2.36 


-3.08 


-3.01 


1.84 


1.79 


-1.18 


-3.38 


-3.37 


1.84 


1.79 


0.00 


-3.73^'') 


-3.73^'') 


1.84 


1.79 


1.18 


-4.07 


-4.10 


1.84 


1.79 


2.36 


-4.42 


-4.46 


1.84 


1.79 


3.53 


-4.77 


-4.82 


1.84 


1.79 


n = no + 2 


7.07 


-1.91 


-2.02 


1.82 


1.79 


8.24 


-2.29 


-2.39 


1.82 


1.79 


9.42 


-2.61 


-2.75 


1.82 


1.79 


10.60 


-2.97 


-3.11 


1.82 


1.79 


11.78 


-3.32 


-3.47 


1.82 


1.79 



('^^ The numbers to be 
electron affinity: -3.31 



compared with experimental values of 
eV Ref. [H] and -3.57 eV Ref. I42j. 



last number is in a very reasonable correspondence with 
the size of the acceptor group of the molecule — see Fig. 

m 

The second and third columns of the table present 
the genuine electron addition energies AE{n) calculated 
from, respectively, Eq. ( A9 1 and (All ), using the average 
of the above values of e^/2C. Not only do these values 
coincide very well; they are in a remarkable agreement 
with experimentally measured electron affinities [HI 
of molecules similar to our molecular trap. 

We believe that these results show that, first, LSDA 
SIESTA provides very small compensation of the self- 
interaction effects in the key energy AE{n) and, second, 
that (at least for the lowest negative ions of our trap 
molecules), an effective compensation may be provided 
using any of the simple relations (A9) and (All I. 



FIG. 18. The schematic energy spectrum of our model at a 
voltage V below voltage Vt that aligns the group localized 
level e with the valence band edge e„. 



occupied level in the group by energy —e{V — Vt), where 
V is the fraction of the voltage drop between the centers 
of the group and the tail of a molecule, and Vt is its value 
which aligns the level with Now let V be close to Vt, 
so that the occupancy p of the discrete level is noticeable. 
If the effect of group charging on the exchange-correlation 
energy is negligible, a generic DFT theory (e.g., LSDA) 
would describe the system energy as 



E = E^-e{V -Vt)v 



(fir 



(fir 



,p(r)p(r') - po(r)po(r') 



(Bl) 



where index marks the variable values bX p — 0. 



Now let us simplify Eq. (Bl ) by assuming that due to a 



small size of the acceptor group, the Coulomb interaction 
of electrons localized on it is much larger than that on 
the chain, so that the latter may be neglected. (For the 
trap molecule shown in Fig. 2, this assumption is true 
within ~5%.) Then Eq. (Im]) is reduced to 



E^Eo-e{V -Vt)p~l f c^{v)\^{v)\^ (fir, 

^ J group 

mvffir, (B2) 



Appendix B: Level freezing in DFT 



group 

where 0(r) is the electrostatic potential created by the 
part of the electronic wavefunction that resides on the 
group. In the simple capacitive model of the group (used 
in particular in Appendix A), 0(r) = —ep/C, where C is 
the effective capacitance of the group, so that 



For the analysis of the fictitious "level freezing" pre- 
dicted by a naive DFT at V > Vt (see Fig. 6|, let us 
consider the following simple model: a molecule consist- 
ing of a small acceptor group with just one essential en- 
ergy level, and a spatially separated chain with a quasi- 
continuous valence band. Figure [18] shows the energy 
spectrum of the system a.t V < Vt- (As before, the occu- 
pied levels are shown in black, while the unoccupied ones 
are shown in green.) 

The edge of the band is separated from the first un- 



E^ Eo-e{V -Vt)p + 



1 2 

e p 



(B3) 



On the other hand, in accordance with Eq. (A7|, the 
total energy in the DFT may also be presented in the 
form 



E = Y^p^e.-^-^'^\(fir\(fir'P^'^P^''^ 



(B4) 
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where Si are all occupied (or partially occupied) single 
particle energies, so that in our simple model 



E Ea + {e - e^)p - 



2 2 

e p 



(B5) 



Comparing Eqs. ( B3 1 and ( B5 ) , we arrive at the following 
expression: 



-e{V - Vt) + 



2 

e p 
~C' 



(B6) 



In most DFT packages, level occupancies pi are calcu- 
lated from the single-particle Fermi distribution, 



(a) 





FIG. 20. (a) A sketch of the evolution of the energy spectrum 
from Fig. [18] as a function of the applied voltage V, illustrat- 
ing the self-interaction errors giving rise to a spurious level 
freezing in the Vt <V <Vt + e/C voltage range. The dashed 
black-green line indicates the region with a partial occupancy 
< p < 1 of the group-localized orbital, (b) A sketch of the 
evolution of the same energy spectrum in a correct quantum- 
mechanical theory, in which electrons do not self-interact. 



exp{(e, -^)/fcBr}-f 1' 



(B7) 



for our simple model, index i may be dropped, and (due 
to the valence band multiplicity) fi Sy. 

As is evident from the sketch of Eqs. (B6) and (B7|, 
in Fig. [19) if the thermal fluctuation scale fceT is much 
lower than the charging energy scale e^/C, then almost 
within the whole range Vt < V < Vt + e/C, the approxi- 
mate solution of the system of these equations is 



p« -(V-I4), e 
e 



(B8) 





P 






1 






Eq. 


(B.6) 









-e(V-V,) n 

FIG. 19. A sketch of Eqs. (IBGIi and (|B7 



Panel (a) in Fig. 20 shows (schematically) the resulting 
dependence of the energy spectrum of our model system 
on the applied voltage V, with level freezing in the range 
Vt < V < Vt+e/C. The dashed black-green line indicates 
the region with a partial occupancy < p < 1 of the 
group-localized orbital. In panel (b) in (Fig. 20 1 we 
show the evolution which should follow from the correct 
quantum-mechanical theory, in which electrons do not 
self-interact and as a result there is the usual anticrossing 



of energy levels e and at V = Vt- (For clarity. Fig. 20 
strongly exaggerates the anticrossing width, which is less 
than 10~^ eV for our trap [5T|.) 



The actual spectrum of our molecular trap is somewhat 
more complex than that of the simple model above — 
see Figs. and |6] First, not only the valence energy 
band of the alkane chain, but also its conduction band 
is important for electron transfer in our voltage range. 
Second, the molecular group has not one, but a series of 
discrete energy levels, with the most important of them 
corresponding to the working orbital (energy £w)i £^nd 
one more group-localized orbital with energy ew+i ~ 
ew + 0-7 eV. 




FIG. 21. (a) A sketch of the evolution of the molecular energy 
spectrum of our trap molecule as a function of the applied 
voltage V, illustrating the self-interaction errors giving rise to 
a spurious level freezing in voltage ranges V/ — e/C <V<Vt 
and Vt < V < Vt+e/C. The dashed black-green line indicates 
the region with a partial occupancy < p < 1 of the group- 
localized orbitals with energies ew or ew+i- (b) A sketch 
of the evolution of the molecular energy spectrum but in a 
correct quantum-mechanical theory, in which electrons do not 
self-interact. 

Nevertheless, the behavior of the spectrum, predicted 
by uncorrected versions of DFT (Fig. [6| may still be 
well understood using our model. Just as was discussed 
above, for voltages V above the threshold value Vt (which 
now corresponds to the alignment of with evv-i-i rather 
than ew), it describes a gradual transfer of an electron 
between the top level of the valence band and the sec- 
ond group-localized orbital, with its occupation number 
Pw+i gradually growing in accordance with Eq. (B8) 
— see panel (a) in Fig. 21 Similarly, at voltages V 
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below VI (which corresponds to the aUgnment of the 
working orbital's energy sw with the lowest level Sc of 
the chain's conduction band), there is a similar spurious 
gradual transfer of an electron between the correspond- 
ing orbitals. In both voltage ranges, a spurious internal 
electrostatic potential is created; as is described by Eq. 
( B8 1 , it closely compensates the changes of the applied 
external potential, thus "freezing" all orbital energies of 
the system at their levels reached at thresholds V/ and Vj 
— see panel (a) in Fig. 21 Figure 22 shows that results 



of both the LSDA and ASIC DFT calculations atV >Vt 
agree well with Eq. (B8), with a value C — 4.5 x 10~^° 



F calculated as discussed in Appendix A, indicating that 
the electron self-interaction effects remain almost uncom- 
pensated in these software packages, at least for complex 
molecules such as our trap. 



Again, in the correct quantum-mechanical theory, 
there should be a simple (and in our molecules, extremely 
narrow) anticrossing between the effective single-particle 
levels of the acceptor group and the alkane chain — see 
panel (b) in Fig. 21 As described in Sec. Ill of the 
main text, we have succeeded to describe this behavior 
rather closely, using the internal iteration dynamics of 
ASIC SIESTA with T = K. 




-Pw?>o) 

-ptv'.T("o + 1) 
-p{^!;?(no + 1) 
liV- V,) 



0.5 1 

V-Vt (Volts) 



1.5 



FIG. 22. The DFT-calculated occupancy pw+i of the {W + 
l)'st orbital of the acceptor group of our trap molecule at 
voltages above the threshold voltage Vt of the alignment of 
its energy ew+i with alkane chain's valence band edge e„. 
Black lines show results of two versions of DFT theory, for two 
ion states: the singly-negatively charged ion and the neutral 
mole cule, while the red line shows the result given by Eq. 



(B8l with C = 4.5 X 10" 
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